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Abstract 

\ The equation of state of the hypothetical dark energy component, which constitutes about two 

\ thirds of the critical density of the universe, may be very different from that of a cosmological 

' constant. Employing a phenomenological model, we investigate semi- analytically the constraints 

^ • imposed on the scalar quintessence by supernovae observations, and by the acoustic scale extracted 

. from recent CMB data. We show that a universe with a quintessence-dominated phase in the dark 

^ ! age is consistent with the current observational constraints. 
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I. INTRODUCTION 



Recent astrophysical and cosmological observations such as dynamical mass, Type 
la supernovae (SNe), gravitational lensing, and cosmic microwave background (CMB) 
anisotropies, concordantly prevail a spatially flat universe containing a mixture of mat- 
ter and a dominant smooth component, which provides a repulsive force to accelerate the 
cosmic expansion The simplest candidate for this invisible component carrying a suffi- 
ciently large negative pressure is a true cosmological constant. The current data, however, 
are consistent with a somewhat broader diversity of such a repulsive "dark energy" as long 
as its equation of state approaches that of the cosmological constant at recent epoch. A 
dynamically evolving scalar field called "quintessence" (Q) is probably the most popu- 
lar scenario so far to accommodate the dark energy component. It is very interesting and 
fundamentally important to distinguish the Q field from the true cosmological constant case. 

Many efforts have been put forth to reconstruct the scalar potential V{(f)) from observa- 
tional data based on various reasonable physical motivations. They include pseudo Nambu- 
Goldstone boson, inverse power law, exponential, tracking characteristics, oscillating feature, 
and others 0]. Several attempts have been made to test different Q-models 3]. Nevertheless, 
it proves to be premature at this stage to perform a meaningful data fitting to a particu- 
lar quintessence model, or to differentiate between the variations. Reconstruction of V{(j)) 
would likely require next-generation observations. 

Since the scalar potential of the Q-field is scarcely known, it is convenient to discuss 
the evolution of (p through its equation of state, = w^p^. Physically, —l<w^< 1, 
where the former equality holds for a pure vacuum state. Lately some progress has been 
made in constraining the behavior of Q field from observational data. A combined large 
scale structure, SNe, and CMB analysis has set an upper limit on Q models with a constant 
W(j) < —0.7 and a more recent analysis of CMB observations gives = — 0.82to;ii 0. 

Furthermore, the SNe data and measurements of the position of the acoustic peaks in the 
CMB anisotropy spectrum have been used to put a constraint on the present < —0.96 0|. 
The apparent brightness of the farthest SN observed to date, SN1997ff at redshift z ~ 1.7, 
is consistent with that expected in the decelerating phase of the fiat ACDM model with 

~ 0.7 Q, inferring w,^ = —1 for z < 1.7. 

Given the above observational constraints, one would like to know the possible role played 
by the Q field in the early universe. If the dark energy is a pure cosmological constant or 
it has a constant equation of state of ^ —0.7, it would be quickly dominated by the matter 
for z > 0.3. The situation may be very different for quintessence with a time-dependent w^. 
Here we will adopt a model-independent approach in which a phenomenological form for 
the time-dependent W(j, is assumed and then used to unfold the dynamics of the Q field up 
to the epoch of the last scattering surface. In particular, we focus on the CMB constraints 
that apply to such a generic quintessence (GQ) scenario in the hope of distinguishing the Q 
field from the true cosmological constant and the constant equation of state. 

II. QUINTESSENTIAL COSMOLOGY 

Consider a fiat universe in which the total density parameter of the universe today is 
represented by fio = + + = 1 with a negligible fi^ and fij^/fij, ~ 1/2. Assuming 
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a spatially homogeneous field, the evolution of the cosmic background is governed by 

3aH{l + w^) g^, (1) 

3 Qj 

-aTi^ - -{WrQr + W^Q^), (2) 

where the conformal time t] = Hq J dta^^it) with the scale factor a and the Hubble constant 
Hq, we have used {dcp/dtY = (1 + w^)pfj, and V^(0) = (1 — w^)p^/2, and rescaled the energy 
density of the i-th component by the reduced Planck mass Mp = (SvrG)"^/^ and Hq, i.e. 
Qi = pi/ [MpHof' . Accordingly, the dimensionless Hubble parameter is given by 

^ (^)' = a-' + a-' + H\ (3) 

Therefore, given a prescribed function of time for the equation of state, w^{r]), the problem 
is reduced to solving a set of first-order coupled ordinary differential equations (P) and 
with the initial conditions set at the present time and the time variable running backwards. 

It is well-known that 7^ after last scattering would result in a time- varying Newtonian 
potential which produces anisotropics at large angular scales through the integrated Sachs- 
Wolfe effect. If we are to constrain dynamical Q models, it will be useful to define an 
fi^-weighted average j9| 

{w^) = n^{T])w^{r])dT] X [ n^{r])dr]\ , (4) 

where r]o and ?2dec are respectively the conformal time today and at the last scattering. It 
was shown |^ that as far as the CMB anisotropy power spectrum is concerned, the time- 
averaged (w^) is well approximated by an effective constant w^p for models in which the Q 
component is negligible at the last scattering. 
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III. THE GENERIC QUINTESSENCE SCENARIO 

The generic quintessence scenario assumes a phenomenological form for the equation 
of state W(j, to accommodate as many observational outcomes as possible. Here we take 
= 0.3 and n° = 9.89 x 10"^ with Wr = 1/3. Note that fi^ = 1 - fi^^ - 1]° ~ 0.7. 
Let us start with the form for at low redshifts. In order to satisfy the above-mentioned 
SNe observational constraints on w<^, we are bound to choose = —1 for z ^ 2. As such, 
the quintessence behaves just like a cosmological constant for < z < 2, and it becomes 
dominated by the matter density at 2; ~ 0.3. This can be seen in an example of the GQ 
in Fig. lU We have also plotted dO/drj = a^J (1 + w^)g^, where 9 = 4>/Mp. For the Q 
field to obtain maximal dynamics at low redshifts, has to increase abruptly to unity 
for z > 2. However, if remains at unity for too long, it would induce an unacceptably 
large integrated Sachs- Wolfe effect on the CMB large-scale anisotropy. In Fig. we have 
chosen a width of the square- wave form of such that {W(j,) ~ —0.7, which saturates the 
above-mentioned upper limit from CMB data on a constant [4]. We thus see that 
can reach a value of 0.42 at z ~ 23 at which is forced to decrease to a negative value. 
On the other hand, a Q component with a constant = —0.7 is dominated by matter for 
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FIG. 1: An example of GQ scenario 

z > 0.5. Also note that at the present time Hot = 0.95, and we have plotted the history of 
the deacceleration parameter q. 

The above GQ model exemplifies an extreme case on how the quintessence may be differ- 
ent from a cosmological constant for z > 2. If we shift the square- wave of to higher red- 
shifts while keeping (w^) ~ —0.7, the integrated Sachs- Wolfe effect on the CMB anisotropy 
power spectrum would be almost the same. Meanwhile, the dynamics of the Q field will 
take place at higher redshifts accordingly. This will be valid until the Q component becomes 
significant at the last scattering. 

IV. CMB ACOUSTIC PEAKS AND QUINTESSENCE 

We now turn to study the effect of the presence of the Q component at the time of 
last scattering on the CMB acoustic peaks. The tightly coupled baryon-photon plasma 
experienced a serial acoustic oscillation just before the recombination epoch. The acoustic 
scale (the angular momentum scale of the acoustic oscillation) sets the locations of the peaks 
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in the power spectrum of the CMB anisotropies and is characterized by 

where represents the comoving distance to the last scattering surface and hs denotes 
the sound horizon at the decouphng epoch with the sound speed Cg. Both d^: and hg are 
affected in the presence of the Q component. Furthermore, the location of the m-th peak 
can be parametrized in practice by the empirical fitting formula, Im = — V'm), where the 
phase shift (frn caused by the plasma driving effect is solely determined by pre-recombination 
physics [11]. It was shown 3] that the shift of the third peak is relatively insensitive to 



cosmological parameters. Consequently, by assuming = 0.341, the value of the acoustic 
scale derived from the analysis of BOOMERANG data 3] hes in the range 3 



/a = 316 ±8, (6) 

which is estimated to be within one percent error if the location ^3 is measured. 

To impose the constraint from the peak separation of CMB, one needs to calculate the 
acoustic scale Ia in Eq. (jSJ for different GO models. Since the sound speed in the pre- 
recombination plasma is characterized by |10| 

^=^= with 



3{1 + R) 



u^^^„,om(^Y^, (7) 

where Hq = 100/ikms^^Mpc^^ and the baryon-photon momentum density ratio R sets the 
baryon loading to the acoustic oscillation of CMB, the sound horizon hg at the decoupling 
epoch can be determined by the differential equation 

dhs 
-d^ = 

coupled with the background evolution equations ((H) and (j2|). Using = 0.05, h = 0.65, 
and fixing the last scattering surface at z^ec = 1100, we have calculated the acoustic scales 
of different GQ models by shifting the square-wave of towards the last scattering surface 
while keeping {w^p) a constant. Figure El plots the results against the Q energy density 
at decoupling Q^{zdec) with contours (w^) = —0.5, —0.6, —0.7, —0.8, —0.9, and —1. The 
last one is simply the ACDM model. In general, when {w,j)) is fixed, the acoustic scale is 
proportional to the quintessence density at the last scattering surface. On the other hand, for 
a fixed fl^{zdec), the acoustic scale decreases as (w^) increases. The two horizontal hues are 
drawn as the upper- and lower-bounds of the acoustic scale derived from the BOOMERANG 
data (jni). All GQ models lying in the region above the curve of {W(f,) = —0.7 and within 
the BOOMERANG bounds are consistent with the current CMB data. Thus, the model 
depicted in Fig. [T] is disfavored as far as the acoustic peak location is concerned. We have 
presented an extreme model located at the upper right corner of the allowed region in Fig. El 
i.e. {{w^) , n^^'^) = (—0.7, 0.59). The figure shows the detailed background evolution of the 
model, in which the Q component dominates over the matter density for 1380 > z > 130 
and the deacceleration parameter q is negative for 1400 > z > 430. 
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FIG. 2: Acoustic scales of the GQ models are plotted as a function of the quintessence density at 
decoupling while keeping the Q^-weighted averaged equation of state fixed. The model equivalent to 
the cosmological constant case has Ia ~ 310 and is denoted by a solid triangle. The two horizontal 
lines signify the upper- and lower-bounds permitted by the BOOMERANG data. 

V. THE AMPLITUDE OF MATTER POWER SPECTRUM 

It is well known that a non-zero would lead to a suppression of the growth of matter 
perturbations. It is therefore useful to study the constraint from the growth of structure 
on the GQ model. Let us consider the growth function for matter perturbations on scales 
smaller than the smoothing scale of the quintessence governed by the differential equation, 

— + an— - -a^n^n^6 = o. 9 

di]"^ drj 2 

We have evaluated the growth function 6{ri) for the extreme GQ model in Fig.lHl Compared 
to the ACDM model, the growth of 6 from the radiation-matter equality to t] = rjo is 
suppressed by about a factor of 2 (see Fig. Ej). When normalized to the COBE data, the 
resulting CMB anisotropy power spectrum has a significant excess of power at / ~ 100. 
In order to keep the suppression of the growth function stable at about 30% level, we 
have found that {wip) has to be less than about —0.9 in Fig. |21 An accurate determination 
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FIG. 3: Background evolution of an extreme GQ model with ((if^) , ^'^^^) = (—0.7, 0.59). The 
last scattering surface is marked as the vertical line at zjec = 1100. 

would require computing small-scale matter perturbations as well as performing a maximum- 
likelihood fitting to the existing CMB anisotropy data. We have plotted in Fig. E] an extreme 
model located at the upper end of the contour {w^) = —0.9 where fl'^^'^ = 0.42, and the 
growth function of this model in Fig. El In this case, the Q component dominates over 
the matter density for 1100 > z > 590 and the deacceleration parameter q is negative for 
1150 > z> 810. 



VI. CONCLUSION AND DISCUSSION 

We have investigated the evolution of the quintessence allowed by the observational con- 
straints from CMB and SNe, using a semi-analytic method with a simple square-wave func- 
tion for the time-varying equation of state. Although the true equation of state, if there is 
any, may be a complicated function of time, the square-wave should roughly capture the 
generic feature of the evolution of the quintessence. This generic quintessence model is 
sufficient for us to confront the current observational data. Future high-precision data will 
tighten the constraints to this model and we may even need a more sophisticated model to 
parametrize the physics of the Q component. Also, the present method gives more physical 
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FIG. 4: Background evolution of an extreme GQ model with {{w,/,) , f^l*"^) = (—0.9, 0.42). 



insights and is much simpler though less accurate than the numerically intensive maximum 
likelihood analysis of CMB data (see, e.g., Ref. jl5l|). 

Three extreme GQ models have been presented. Figure shows the maximum dynamics 
that the Q field can attain at low redshifts for z > 2. The evolving Q field during the 
large-scale-struture formation may have interesting cosmological implications. For instance, 
the authors in Ref. M have attempted to generate primordial magnetic fields from the 
dynamics of the Q field coupled to electromagnetism. This electromagnetic Q field may also 
be responsible for the time- varying fine structure constant (a) jT3| as it was recently claimed 
that the results of a search for time variability of a using absorption systems in the spectra 
of distant quasars yield a smaller a in the past 0]. We have also studied the constraint 
from the growth of large-scale matter perturbations on the GQ model. Figure |3 shows 
that the Q component can make up about 40% of the total energy density of the universe 
at last scattering. This result is consistent with the upper bound fl^ < 0.39 during the 
radiation dominated epoch obtained by performing a maximum likelihood analysis on the 
CMB data In general, the GQ scenario bears a salient feature that the Q component 
overwhelms the matter during the dark age. It is worth studying in more details about its 
influence on the evolution of matter perturbations and the subsequent structure formation. 
At last, we would like to point out that an acceleration of the universe in the past is consistent 
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FIG. 5: Growth functions of large-scale matter perturbations with unity normalization at radiation- 
matter equality time for the ACDM model and the GQ models in Fig. and Fig. ^ 

with all observations. So the fact that the universe is accelerating today would not be quite 
unnatural. 
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